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THE  MEASUREMENT  OP  AERODYNAMIC  FORCES  AND  MOMENTS 
IN  THE  NOL  4-IN.  HYPERSONIC  SHOCK  TUNNEL  NO.  3 

Prepared  by: 

D.  P.  Gates  and  D.  N.  Bixler 


ABSTRACT:  Two  light-weight  models  of  missiles  are  suspended  by 
very  fine  threads  in  the  test  section  of  the  NOL  4-ln.  Hypersonic 
Shock  Tunnel  No.  3.  Alongside  the  models,  a  light-weight  sphere 
is  suspended.  A  high-speed  camera  is  focused  on  the  models 
throxigh  a  window  in  the  side  of  the  shock  tunnel  receiver  tank 
and  lll\jminatlon  of  the  models  is  provided  by  a  short -duration, 
hlgh-lntensity  light  source.  When  the  shock  tunnel  is  fired, 
the  light  source  is  actuated,  the  fine  threads  holding  the 
models  are  snapped,  and  the  models  are  then  photographed  for 
2-1/4  milliseconds  as  they  move  freely  in  the  shock  tunnel  flow. 
Prom  the  high-speed  camera  film  (taken  at  a  frequency  of  about 
one  frame  every  40  microseconds)  data  of  the  horizontal,  vertical, 
and  angular  motions  of  the  models  as  functions  of  time  are 
obtained.  Data  of  the  motion  of  the  sphere  are  also  obtained. 
Since  the  aerodynamic  characteristics  of  a  sphere  are  fairly 
well  established,  this  configuration  is  employed  as  a  calibration 
body,  and  the  motions  of  the  models  are  compared  to  that  of  the 
sphere.  By  various  methods  of  comparison  the  aerodynamic  drag, 
lift,  and  moment  coefficients  of  the  models  are  determined. 
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THE  MEASUREMENT  OP  AERODYNAMIC  FORCES  AND  MOMENTS 
IN  THE  NOL  4-IN.  HYPERSONIC  SHOCK  TUNNEL  NO.  3 


INTRODUCTION 

1.  In  order  to  create  conditions  In  the  laboratory  similar  to 
the  conditions  encountered  by  a  high-speed  missile  in  a  re-entry 
trajectory  at  an  altitude  of  50,000  to  100,000  feet,  the  Naval 
Ordnance  Laboratory  constructed  the  NOL  4-In.  Hypersonic  Shock 
Tunnel  No.  3-  On  3  September  1958,  the  first  experiment  with 
this  facility  was  performed.  Since  then,  most  of  the  tests  In 
this  shock  tunnel  have  been  made  for  the  purpose  of  determining 
the  static  aerodynamic  drag,  lift,  and  moment  coefficients  of 
models  of  missiles  as  functions  of  missile  angle  of  attack  and 
Mach  number.  In  each  test,  this  was  accomplished  by  observing 
the  motion  of  light-weight  models  in  the  test  section  while 
these  models  were  In  actual  free  flight  In  the  shock  tunnel  flow. 
The  experimental  technique  followed  in  each  test  was  first 
briefly  described  by  representatives  of  the  Naval  Ordnance 
Laboratory  at  the  Second  Shock  Tube  Symposium  In  1958  (reference 
(a)).  A  more  elaborate  description  was  given  in  April  1959 
(reference  (b))  and  in  May  1959  (reference  (c)).  The  technique, 
with  several  modifications.  Is  also  being  used  by  the  Missile 
and  Space  Vehicle  Department  of  the  General  Electric  Company 
In  Its  6-ln,  shock  tunnel  (reference  (d)).  Since  the  first  test 
was  made  In  1958,  different  experimental  procedures  and  methods 
of  data  analysis  have  been  employed.  All  of  the  various  experi¬ 
mental  and  analytical  methods  are  described  below. 


EXPERIMENTAL  TECHNIQUE 

2.  Figure  1  is  an  Illustration  of  the  NOL  4-ln.  Hypersonic 
Shock  Tunnel  No.  3.*  Figures  2  and  3  give  some  of  the  dimensions 
of  the  various  parts  of  the  shock  tunnel.  The  high-pressure 
chamber  (which  subsequently  will  be  referred  to  as  the  "chamber") 
Is  12  feet  long  and  has  an  Inside  diameter  of  10  Inches.  The 
low-pressure  chamber  (which  will  be  referred  to  as  the  "barrel") 
Is  60  feet  long  and  has  an  Inside  diameter  of  4  Inches,  while 
the  receiver  tank  Is  54  feet  long  and  8  feet  In  diameter. 

Between  the  chamber  and  the  barrel  Is  a  flat  stainless-steel 
diaphragm  which  has  a  thickness  of  about  0.3  Inches.  A  piston 
Is  placed  next  to  the  diaphragm,  though  In  many  tests  the 


*  A  general  description  of  this  facility  and  Its  Instrumentation 
can  be  found  In  reference  (e). 
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piston  Is  excluded.  It  Is  made  of  nylon  or  aluminum  and  weighs 
from  200  to  550  grams.  Several  pressure  gages*  are  mounted 
through  holes  at  various  stations  along  the  side  of  the  barrel, 
so  that  the  pressure  Pm  of  the  air  at  these  stations  can  be 
monitored.  At  the  receiver  tank  the  barrel  diameter  Increases 
conically  from  4  inches  to  10  Inches.  This  Increase  takes 
place  over  a  distance  of  9.6  Inches,  and  at  a  distance  of  6.4 
Inches  from  the  large  end  of  the  Increase  there  Is  a  hemispherical 
steel  diaphragm  6  Inches  In  diameter  and  tapered  In  thickness 
from  0.12  Inch  at  the  center  to  0.30  Inch  at  the  perimeter. 

This  muzzle  arrangement  serves  to  prevent  disturbances  from  the 
opening  hemispherical  diaphragm  from  Interfering  with  the  flow 
over  the  models  by  allowing  these  disturbances  to  diverge  as 
they  enter  the  test  section. 

3.  The  flow  Is  allowed  to  expand  In  a  free  Jet  from  the  muzzle 
Into  the  receiver  tank.  Thus,  the  flow  over  the  models  In  the 
test  section  Is  radial,  and  the  test  section  Mach  number  Increases 
with  Increasing  distance  from  the  muzzle.  At  various  distances 
from  the  muzzle,  pairs  of  opposing  windows  are  located  In  the 
sides  of  the  receiver  tank,  as  shown  In  Figure  3>  so  that  models 
may  be  placed  between  any  pair  of  windows,  depending  upon  the 
flow  conditions  desired.**  The  distances  between  the  muzzle  and 
the  windows  shown  In  Figure  3  may  be  conveniently  changed  by 
adding  a  3-foot  extension  to  the  barrel.  Also,  the  most  distant 
three  windows  are  all  mounted  on  a  plate  which  can  be  shifted 
downstream  by  one  foot.  Models  are  most  commonly  tested  at 
distances  to  80,  128,  or  2l8  Inches  from  the  muzzle.  Table  I 
lists  the  standard  loading  conditions  and  approximate  flow 
conditions  observed  when  testing  at  these  positions  In  the  test 
section.  The  flow  duration  Is  approximately  2-1/4  milliseconds, 
except  for  the  Mach  l4.3  case,  where  It  Is  slightly  less  than 
2  milliseconds.  In  all  cases  shown  In  this  table,  the  chamber 
Is  loaded  with  6,000  psl  of  helium,  and  the  hellimi  Is  heated  by 
burning  30  lbs.  of  Hercules  Unique  powder.  The  powder  Is  Ignited 


*  A  detailed  description  of  the  type  of  gages  used  In  this  shock 
txmnel  Is  given  In  references  (f)  and  (g). 

♦♦Tests  have  shown  that  between  any  pair  of  windows  the  axial 
gradient  of  the  Mach  number  over  a  distance  of  4  inches  Is 
always  less  than  2-1/2  percent  of  the  Mach  number,  while  cal¬ 
culations  have  revealed  that  radial  gradients  In  Mach  nvunber 
over  a  distance  of  6  Inches  off  the  centerline  of  flow  are 
always  less  than  2-1/2  percent  of  the  Mach  number  (cf., 
reference  (h)). 
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by  a  series  of  explosive  primers  running  the  length  of  the 
chamber.  In  some  Instances  the  chamber  is  loaded  with  a  mixture 
f  of  oxygen,  hydrogen,  and  helium,  in  a  l/3/B  molar  proportion, 

respectively.  Instead  of  helium  and  powder.  This  mixture  is 
also  ignited  by  the  same  arrangement  of  explosive  primers. 

Though  the  helium  powder  mixture  produces  a  less  efficient 
driver,  it  is  more  commonly  used  because  there  is  little  chance 
of  an  undesirable  detonation  occurring  during  ignition  of  the^ 
powder.  At  the  same  time  efficiency  is  not  the  object,  but 
high  densltltes  and  long  flow  durations  are  more  desirable. 

4.  In  the  test  section,  two  light-weight  models  and  a  light¬ 
weight  sphere  are  suspended  by  very  fine  threads.  Figure  4  is 
a  photograph  of  two  models,  (in  this  case,  hollow  cones)  and 

a  sphere  suspended  by  threads.  This  is  a  view  from  the  rear  of 
the  receiver  tank  looking  forward  into  the  muzzle.  Outside  one 
window  a  short -duration,  high-intensity  light  source  is  situated 
so  as  to  pass  a  converging  beam  of  light  through  the  receiver 
tank  and  over  the  models.*  A  system  of  several  horizontal  and 
vertical  grid  wires  is  attached  to  the  window  nearest  the  light 
V  source.  This  system  of  wires  serves  as  a  reference  system,  and 

can  be  removed  from  one  window  and  placed  on  any  other  window. 

A  Beckman  &  VHiltley  Dynafax  high-speed  framing  camera  is  placed 
on  the  opposite  side  of  the  receiver  tank  and  focused  on  the 
models.  This  camera  is  capable  of  taking  pictures  at  a  rate 
as  high  as  one  frame  per  40  microseconds.  In  tests  where  a 
greater  framing  rate  is  desired,  the  Beckman  &  Whitley  192  model 
high-speed  framing  camera  is  used;  however,  inherent  in  this 
camera  is  a  time-varying  optical  distortion  (cf.,  reference  (j)). 
The  Dynafax  model  also  has  a  minor  amount  of  distortion,  but  it 
is  essentially  constant  with  time. 

5.  The  models  suspended  in  the  test  section  weigh  between  one 
and  30  grams,  and  are  from  one  to  6  Inches  long.  They  can  be 
manufactured  either  by  machining  or  by  molding.  If  they  are 
machined,  they  are  usually  made  of  ethocel,  althougn  various 
materials  are  sometimes  used  for  different  parts  of  a  model. 

A  portion  may  be  made  of  magnesium,  aluminum,  brass,  or  steel, 
in  order  to  Increase  strength  of  a  particular  part  of  a  model 
or  to  position  the  center  of  gravity  at  a  desired  point.  The 
tolerance  on  the  machining  is  either  +  0.002  inch  or  +  O.OO5 
inch,  depending  upon  the  material  macKlned.  When  a  model  is 


*  A  detailed  description  of  the  light  source  is  given  in 
reference  (l). 
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of  such  a  design  that  It  cannot  be  machined  Inexpensively,  while 
at  the  same  time  a  sufficiently  large  quantity  of  models  Is 
required,  a  mold  Is  built,  and  the  models  are  then  molded.  In 
this  case,  they  are  manufactured  from  a  foamed  plastic,  and  the 
tolerance  on  each  model  Is  kept  within  +  O'. 002  Inch.  The  sphere 
that  accompanies  the  models  In  the  test  section  has  a  diameter 
of  0.25  to  1.5  Inches,  the  latter  being  the  diameter  of  a  ping 
pong  ball,  which  Is  sometimes  used.  Each  sphere  weighs  from 
0.3  to  5  grams  and,  except  for  the  ping  pong  ball.  Is  made  of 
nylon. 

6.  With  two  models  and  a  sphere  suspended  In  the  test  section 
and  the  shock  tunnel  Instrvimentatlon  prepared,  the  test  section 
Is  sealed  off  and  evacuated  to  less  than  30  microns  pressure. 

The  chamber  Is  loaded  with  30  lbs.  of  powder,  the  primer  rod 
Is  Inserted,  and  the  breech  Is  secured.  The  barrel  Is  loaded 
with  the  desired  pressure  of  air,  while  the  chamber  Is  loaded 
with  helium.  The  camera  Is  then  started,  and  when  the  camera 
drum  attains  a  rate  of  rotation  of  about  1,560  rps,  the  camera 
shutter  Is  opened,  and  the  primers  are  electrically  detonated. 
Igniting  the  powder.  The  hellvim,  due  to  the  powder  combustion. 

Is  raised  to  a  high  temperature  and  pressure,  causing  the 
stainless-steel  diaphragm  to  rupture.  The  heated  hellm  then 
drives  against  the  air  In  the  barrel,  forming  a  shockwave  which 
propagates  through  the  air  toward  the  muzzle,  raising  the 
temperature,  pressure,  and  velocity  of  the  air,  until  the  shock- 
wave  strikes  the  hemispherical  diaphragm,  bursting  It  and  allowing 
the  air  to  expand  Into  the  test  section  and  pass  over  the  models. 
The  light  source  is  actuated  by  a  piezoelectric  shockwave  detector* 
In  the  muzzle,  and  the  threads  holding  the  models  are  swept  away 
by  the  Initial  portion  of  the  flow.  The  models.  Initially  at 
rest,  are  then  accelerated,  and  as  they  move,  they  are  photographed 
by  the  high-speed  camera.  Prom  the  camera  film,  data  of  the 
movements  of  each  model  are  recorded,  and  the  aerodynamic  coef¬ 
ficients  of  drag,  lift,  and  moment  of  the  models  are  reduced  from 
the  data. 


DATA  ANALYSIS  • 

7.  All  methods  of  data  analysis  developed  herein  are  derived 
on  the  assumption  that  dynamic  forces  and  moments,  as  well  as 


*  The  shockwave  detector  Is  described  In  reference  (k). 
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the  force  of  gravity*,  are  negligible.  Thus,  it  is  assumed 
throughout  that  the  magnitudes  of  the  various  components  of 
acceleration  of  a  model  are  directly  proportional  to  the  static 
aerodynamic  moment,  lift,  and  drag  on  the  model.  The  drag  force 
on  a  model  is 


(1) 


where  Is  the  model's  drag  coefficient,  Qm  Is  Its  dynamic 
pressure,  Its  characteristic  area,  m  Its  mass,  and  x  Its 
acceleration  In  the  direction  of  the  flow.  The  lift  on  the 
model  Is 


C,q«Am=my 


(2) 


where  Gl  Is  the  lift  coefficient  of  the  model,  and  y  Is  its 
acceleration  normal  to  the  direction  of  flow  and  In  a  plane 
parallel  to  the  plane  of  the  camera  film.  The  static  aerod^Tnamlc 
moment  on  the  model  Is 


qm  Am  Dm  I  6<  (  3  ) 


Where  Cm  Is  the  static  moment  coefficient  of  the  model,  Dj^  is 
the  characteristic  diameter  of  the  model,  I  is  the  moment  of 
Inertia  about  Its  center  of  gravity,  and  6(  is  its  angular 
acceleration.  Similarly,  for  the  sphere  in  the  test  section 
the  drag  Is 


C,^q,As="msX.s  (^) 

where  Cpg  Is  the  drag  coefficient  of  the  sphere,  pg  is  the 
dynamic  pressure  on  the  sphere,  Ag  is  its  cross-sectional  area, 
mg  Its  mass,  and  Xg  Its  acceleration. 


♦  The  force  of  gravity  on  a  model  is  so  small  In  comparison  to 
the  aerodynamic  forces  that  its  effects  are  undetectable  in 
the  short  dviratlon  in  which  observations  are  made. 
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8.  By  various  manipulations  of  equations  (l)  through  (4), 
values  of  Cp,  Ct^  and  Cjvj  are  determined  directly  from  the 
experimental  data.  Through  equation  (4)  the  dynamic  pressure 
experienced  by  each  model  is  determined.  Since  the  loading  - 
pressure  and  temperature  and  final  pressure  Pf/j  in  the  muzzle 
are  known  from  direct  measurements,  the  entropy  of  the  air  can 
be  found  by  use  of  the  equations  of  conservation  of  fluid  mass, 
momentum,  and  energy,  and  the  equation  of  state  of  air.*  It  is 
assumed  that  the  expansion  of  the  air  from  the  muzzle  is  isentropic, 
so  that  a  value  for  the  entropy  at  any  point  in  the  test  section 

is  known.  Then  the  Mach  number  and  various  other  properties  of 
the  flow  in  the  test  section  can  be  calculated  by  use  of  the 
same  equations  of  fluid  dynamics  and  the  equation  of  state,  once 
the  dynamic  pressure  q  is  found.  For  this  purpose,  equation  (4) 
is  again  used  in  conjunction  with  the  experimental  data  to 
determine  the  dynamic  pressure.  In  these  applications  of  equation 
(4),  two  assumptions  are  made:  (1)  the  dynamic  press\ire  qs  on 
the  sphere  is  equal  to  the  dynamic  press\ire  qjj,  on  the  model,  and 
(2)  the  dynamic  pressure  qg  on  the  sphere  is  equal  to  the  dynamic 
pressure  q  as  experienced  By  an  observer  standing  in  the  test 
section. 

9.  After  a  finite  period  of  time  from  the  beginning  of  flow, 
however,  the  models  and  the  sphere  will  have  attained  individual 
velocities,  so  that  each  will  experience  a  dynamic  pressure 
slightly  different  from  that  experienced  by  the  others,  and  all 
will  experience  a  dynamic  pressure  different  from  that  experienced 
by  a  stationary  observer  in  the  test  section.  To  an  observer 
standing  in  the  test  section,  the  dynamic  pressure  is  q=l/2y9Ll^, 
where  p  and  u  are,  respectively,  the  density  and  the  velocity 

of  the  flow.  Prom  the  standpoint  of  the  sphere,  q3=l/2  p  (u-Ug)^, 
where  Ujls  the  velocity  of  the  sphere.  On  the  model, 

(u-Ujjj)2,  where  um  is  the  velocity  of  the  model.  Now,  tne  dif¬ 
ference  in  dynamic  pressures  as  seen  from  the  model  and  as  seen 
from  the  sphere  is 

q»-q.  =y/>(u-u„;'-i/j(u-u.y 

=^/>[(U-uJ-(u-Us/] 

=5/3{(u-u„)'-  [(u„-  u,) +(u  - 

=iyo[2(u.-uJ(U-a.)+(u„-Usy]  (5) 

•Values  for  the  equation  of  state  of  air  are  taken  from  references 
(1)  and  (m). 
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The  difference  In  dynamic  pressures  as  seen  from  a  point 
stationary  in  the  test  section  and  as  seen  from  the  sphere  is 

p(u-u.y 

=i/ofu'-(u-Us/] 

Usj  (5^ 


The  models  or  the  sphere  seldom  reach  a  velocity  greater  than  2 
percent  of  the  flow  velocity.  In  the  most  extreme  case,  the 
sphere  attained  a  velocity  Ug  of  300  ft/sec,  but  since  the 
velocity  of  the  flow  is  6,300  ft/sec  or  higher,  Ug  and  are 
still  less  than  5  percent  of  u.  Thus,  equations  (5)  and  (6)  can 
be  simplified  by  dropping  out  certain  second-degree  terms  in 
Um  and  Ug.  These  equations  then  become 

q«  -  q,  ^  ~p(u„-  Us) (u  -  U„,) 

q-qs="yOuUs 

The  percent  errors  introduced  In  the  quantities  Qm  and  q  by 
observing  the  dynamic  pressure  from  the  standpoint  of  the  sphere 
are  then 


S'"  ^yO(U-U„) 

=  (7) 

100 -4^' 100 

(8) 

In  each  Individual  test,  the  models  and  the  sphere  are  so  designed 
that  the  difference  between  the  velocity  of  each  model  and  that  of 
the  sphere  (l.e.,  u^  -  Ug)  Is  less  than  one  percent  of  u  -  10^. 
Thus,  according  to  equation  (7),  the  maxlm\jm  difference  between 
qg  and  q^  is  less  than  2  percent  of  q^.  Therefore,  the  first 
assumption,  qm  =  qs,  involves  a  negligible  error.  The  second 
assumption  is  that  q  =  qg.  Thus,  according  to  equation  (8), 
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the  error  resulting  from  this  assumption  Is  less  than  10  percent 
of  q,  since  Ug  Is  always  less  than  5  percent  of  u.  This  . 

assxAmptlon  Is  made  for  the  purpose  of  conveniently  determining  * 

the  Mach  n\;mber  of  the  flow,  and,  as  pointed  out  In  reference 
(b),  a  10  percent  error  In  determining  the  dynamic  pressure  q 
results  In  an  error  of  less  than  2  percent  In  Mach  number. 

With  the  use  of  the  above  assiamptlons,  the  dynamic  pressure  can 
now  be  determined  by  analysis  of  the  sphere  motion  as  photographed 
by  the  camera. 

10.  Obtained  from  the  camera  Is  a  strip  of  film  containing  30 
to  80  consecutive  frames  of  silhouettes  of  the  models  and  the 
sphere  In  motion.  This  film  Is  placed  In  the  projector  of  the 
Universal  Telereadex,*  and  Images  of  the  silhouettes  In  one 
frame  are  magnified  and  cast  on  an  opaque  screen.  A  full-scale 
drawing  or  template  of  the  image  of  each  model  Is  made  on  a  sheet 
of  vellum,  with  the  centerline  and  center  of  gravity  of  the  model 
also  placed  on  each  template.  These  templates  are  then  super¬ 
imposed  upon  the  Images  of  the  silhouettes  In  each  consecutive 
frame  of  the  film,  and  with  the  read-out  system  of  the  Telereadex, 

the  coordinates  of  the  center  of  gravity  and  one  other  position  * 

on  the  centerline  of  each  model  are  recorded  on  IBM  cards.  Also 
recorded  In  the  same  manner  are  the  coordinates  of  the  sphere  In 
each  frame.  As  was  mentioned  earlier,  however,  the  images  In  « 

each  frame  of  the  film  are  distorted  to  some  degree  by  the 
optical  system  of  the  high-speed  camera. 

11.  If  the  Beckman  &  Whitley  192  camera  Is  used,  then  the  film 
taken  from  the  camera  Is  35-nim  size.  Due  to  an  Inherent  property 
of  the  optical  system  of  the  camera.  Images  In  each  frame  of  the 
film  are  distorted  differently  from  those  of  any  other  fraune.  To 
compensate  for  this  distortion,  a  map  of  a  fine  grid  system  Ir 
each  frame  of  the  film  Is  superimposed  on  an  xondlstorted  cartesian 
coordinate  system.  Then  when  an  arbitrary  point  In  any  given 


*  Model  29A,  built  by  Telecomputing  Corporation,  consists  of  a 
projector  which  magnifies  film-size  images  as  much  as  44-1/2 
times  and  casts  these  Images  on  a  white  opaque  screen.  On  the 
screen  are  a  pair  of  perpendicular  cross-wires,  which  can  be 
manipulated  by  an  operator.  Magnetic  reading  heads  In  the 
machine  sense  the  movements  of  these  cross-wires  and  convey  this 
Information  to  the  Telecordex,  which  records  and  reads  out  the 
position  of  the  Intersection  of  the  cross-wires.  The  magnetic 
reading  heads  sense  a  movement  of  the  cross-wires  of  1  Inch. 

395 
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frame  Is  observed  In  the  cartesian  coordinate  system  of  the 
Universal  Telereadex,  the  observed  coordinates  are  replaced  by 
the  corresponding  coordinates  In  the  mapped  distorted  system, 
and  the  values  of  the  latter  coordinates  are  taken  as  the  actual 
position  of  the  point. 


12.  If  the  Beckman  &  Whitley  Dynafax  Is  used,  the  models  are 
photographed  on  l6-mm  frames,  situated  alternately  on  each  edge 
of  a  strip  of  35-nim  film.  Images  In  the  frames  on  one  edge  of 
the  film  suffer  a  small  amount  of  distortion,  while  Images  In 
the  frames  on  the  opposite  edge  suffer  a  different  distortion. 
However,  the  distortion  In  the  frames  on  each  edge  of  the  film 
Is  nearly  constant  with  time.  When  data  of  the  position  of  a 
model  with  respect  to  time  are  taken  from  the  film,  the  data 
points  appear  alternately  deflected,  as  a  result  of  the  alternate 
distortion.  In  order  to  reduce  the  effects  of  this  distortion, 
every  coordinate  In  each  frame  of  the  data  Is  averaged  with  the 
corresponding  coordinates  In  the  adjacent  frame.  For  example. 

If  In  frame  number  n  the  coordinates  of  the  center  of  gravity  of 
a  model  are  Xn  and  y^,  and  In  frame  number  n  +  1,  the  coordinates 
are  Xq  4.  i,  and  yn  +  1  then  the  final  value  given  for  the  position 
of  the  model  Is 


0 


The  data  obtained  from  the  film  are  thus  reduced  to  the  graphical 
form 

y=y(0 

o(=o<(0  ,  ,  ^ 

Xs-Xs(0  ^ 


where  x  and  y  are  the  coordinates  of  the  center  of  gravity  of 
the  model,  x  Is  in  the  direction  of  flow  at  the  centerline  of 
flow,  while  y  Is  normal  to  x  and  In  the  plane  of  the  film. 
Similarly,  Xg  and  ys  are  the  coordinates  of  the  center  of  gravity 
Of  the  sphere.  The  angle  of  attack  oc  of  the  model  Is  found  from 
the  relation 


o  =arc  tan 


where  xg  and  y2  are  the  coordinates  of  any  other  point  on  the 
centerline  of  the  model.  All  coordinates  are  taken  with  respect 
to  some  fixed  arbitrary  reference  point  In  the  grid  system  which 
appears  In  each  frame  of  the  film. 
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13.  Since  the  flow  of  air  in  the  test  section  is  radial,  the 
observed  motion  of  a  model  will  be  slightly  In  error.  This  is 
because  the  x-axls  is  set  parallel  to  the  centerline  of  flow 
and  the  motion  of  the  model  is  then  assumed  to  be  under  the 
Influence  of  a  flow  that  is  everywhere  parallel.  To  correct  this 
discrepancy  the  following  procedure  is  used.  In  earlier  tests 
it  was  estimated  that  a  distance  X  from  the  diaphragm  at  the 
muzzle  along  the  centerline  of  flow  in  the  test  section,  the 
direction  of  the  flow  at  some  point  a  small  distance  Y  from 
the  centerline  is  such  that  the  flow  appears  to  have  originated 
at  a  point  on  the  centerline  a  distance  approximately  0.1  X 
from  the  diaphragm.  Thus,  when  a  model  is  at  a  distance  X  from 
the  diaphragm  and  a  distance  Y  from  the  centerline,  the  observed 
motion  of  the  model  may  be  corrected  by  subtracting  the  quantity 
X  Y/0.9  X  from  the  y-coordinate  of  the  model.  The  quantity  Y  is 
determined  by  observing  the  motion  of  the  sphere.  A  linear 
equation  of  the  form  y^  =  fitted  by  the  method  of 

least  squares  (see  Appendix)  to  the  x^  versus  yg  data,  and  the 
first  derivative  of  this  equation  is  foiind.  On  the  assump- 

tlon  that  the  sphere  travels  in  the  direction  of  the  flow,  the 
distance  Yg  from  the  sphere  to  the  centerline  of  flow  is  then 
taken  as  Then  Y=Yg+y-yg.  As  a  result  of  the  very 


small  Imperfections  in  its  mass  distribution  and  dimensions,  the 
sphere  does  not  travel  precisely  in  the  direction  of  flow. 

However,  the  amount  that  the  sphere  trajectory  deviates  from 
the  direction  of  flow  is  usually  small  compared  to  the  angle 
between  the  flow  direction  and  the  horizontal.  In  one  exception, 
where  x  i  193  Inches  (the  region  in  which  the  flow  Mach  number 
is  greater  than  13-5),  the  magnitude  of  the  flow  radiallty  is 
less  than  the  probable  error  Introduced  when  assuming  that  the 
sphere  travels  in  the  flow  direction.  Therefore,  in  this  region, 
no  correction  is  made  for  flow  radiallty.  For  the  sake  of 
simplicity,  it  will  henceforth  be  assumed  that  the  data  describing 
equations  (9)  have  been  corrected  for  flow  radiallty  in  the 
region  where  this  correction  is  appropriate. 

l4.  The  problem  remaining  is  to  determine  the  quant itltes 
C-Q,  Cl,  and  Cj^^j  as  functions  of  cx  and  Mach  nimiber  M,  by 
manipulating  equations  (l)  through  (4)  and  the  experimental  data 
describing  equations  (9).  Since  the  first  tests  were  made  in 
the  NOL  4-in.  Hypersonic  Shock  Tunnel  No.  3*  various  techniques 
have  been  employed  to  obtain  the  aerodynamic  coefficients.  Each 
method  is  described  in  the  outline  below.  Each  of  Methods  I,  II, 
and  III  is  a  technique  for  determining  Cp,  CL,  and  Cj^,  while  Methods 
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IV  and  V  are  techniques  used  for  determining  values  for 
only.  Methods  VI  and  VII  are  procedures  for  calculating  the 
dynamic  pressure  q. 

15.  Method  I:  By  dividing  each  of  equations  (l)  through  (3)  by 
equation  (4),  and  manipulating  the  results,  one  can  obtain  the 
following  three  relations. 


n  -r 

^0  M,.;  A  fn.  X, 

/  'm  ‘  • 


(10) 


(11) 


~  c''  ^  J  cx 


(12) 


The  quantity  Cpg  Is  a  fairly  well  established  quantity,  and  is 
considered  to  be  equal  to  0.88  in  the  realm  of  Mach  numbers  8 
to  15.  The  other  quantitites  on  the  right-hand  side  of  the 
above  three  equations  are  readily  measurable,  except  the  ratios 
of  accelerations  x  ,  y  ,  and  «  .  These  ratios  are  obtained  by 

^s  ^s  ^s 

first  dividing  the  data  of  equations  (9)  into  relatively  short 
intervals  of  time  and  assuming  that  over  each  time  interval 
accelerations  are  constant.  A  second-degree  equation  is  then 
fitted  to  each  set  of  the  data  by  the  method  of  least  squares 
(see  Appendix),  and  each  equation  is  differentiated  twice  with 
respect  to  time.  This  second  derivative  is  then  taken  as  a 
value  for  the  acceleration.  For  example,  an  equation  of  the 
form  Of  =CQ+Cit+C2t‘^  is  fitted  to  the  0  =  5^  (t)  data,  while  an 


equation 
=  Xs  (t) 


of  the  form  X3=Co, 
data.  These 


+  Cjgt  +  C23t‘^  is  fitted  to  the  X, 


two  equations  are  then  differentiated 


twice,  and  the  ratio  of  these  derivatives. 


the  ratio  of  accelerations*^ 


is  set  equal  to 


Xs 


The  moment  coefficient  is 


then  determined  through  the  equation. 


n  -  C 

A  rn,Dm  C, 
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The  drag  and  lift  coefficients  are  also  determined  In  a  similar 
manner.  Thus,  by  this  method  of  analysis,  all  three  aerodynamic 
coefficients,  C^,  Cl»  and  Cj),  of  a  model,  can  be  determined  If 
the  model  and  tne  sphere  accelerate  In  a  quasi-steady  manner. 
This  condition  Is  met  when  the  angle  of  attack  o<  and  dynamic 
pressure  q  experience  only  small  changes  In  a  given  Interval 
of  time. 

l6.  Method  II:  A  double  integration  of  equation  (l)  yields 


where  x=X4  and  dx  =  0  when  t=0.  is  defined  by 


C- 


JtCpqdtdt 


//qdtdt 

o  o  * 

A  similar  double  Integration  of  equation  (4)  yields 


(14) 


dx. 


where  Xg=x-,4  and  ^s=  q  when  t  =  o.  Cp.  is  assumed  to  be 

SI  ^  Ug 

constant.  Dividing  equation  (13)  by  equation  (l4)  gives 

t  t 

X:2ii=:^!2kS» 

xrXsi  As  rn  Ci 


and  finally 


“  -vv- - 


A, 


p  ^  p  ZlkiZL  X-Xi 
^  Am  Ws  Xs“ Xjsi 


(16) 


Analogously,  if  equations  (2)  and  (3)  are  Integrated  twice, 
and  each  resulting  equation  divided  by  equation  (15),  one  will 
obtain 


CrC. 


A 

A 


s  fH 

rris 


YrJL 

Xs-Xsl 


(17) 


I 
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C.rC, 


I 


0<-0<j 


M  A™  DmlDs  ^S-Xsi 
Where  7=7^,  cx  ^ 


a 

dt 


(18) 

=  0  when  t  =  0,  and  and 


are  defined  b7  the  relations 


^  _  JjQqdtdt 
u,  "=  -y% - 

Jfqdtdt 


(19) 


ffqdtdt 


(20) 


As  cx  and  q  approach  constant  values  in  the  entire  time  interval 
zero  to  t,  the  values  of  and  TTjvi  approach  Cq,  Cl  and  Cjvi, 

respectlvel7. 


17.  Method  III;  If  a  model  is  axiSTinmetric ,  its  moment  and  lift 
coefficients  are  odd  functions  of  its  angle  of  attack  or  .  Thus, 


Cj“  Acx  +  Bcx^  +  Co^  -h  ■  •  • 

Cl^AcX  +  B'o»^  +  C'cx^  +  •  ■ 

At  the  same  time,  the  drag  coefficient  of  an  axisymmetrlc  model 
is  an  even  function  of  cx  »  so  that 

Cd=Co+Co,«’’+  Co^o<’  + 

If  or  is  ver7  small,  then  one  ma7  accuratel7  assume  that 
Cl=A'o(  ,  and  Cq  =  Cq,  .  For  Intermediate  angles  (-50°<cx<  50°) 
reasonable  approximations  are 


Cm=Ao<  +  Bo('^ 


(21) 
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Cf  A'o<  +  Ba’  (22) 

C=C„+  Q„cx=  (23) 


Q 


Upon  substituting  equations  (21),  (22),  and  (23)  Into  equations 
(20),  (19),  and  (l4),  respectively,  one  obtains 


1 1  t  t  t  t 

Cjjqdt  dt = Ajj^qdt  dt  +  bJ jcx^dt  dt  ^  ^ 

o  o  o  o  o  o 


•t  t  t  t  1 1 

Cjjqdtdt  -  A'Jj«c]dt  d  t  +  t  dt  ( 25 ) 

o  o  o  o  o  o 


1 1  t  i  1 1 

CoJJqdtdt  =  CoJJqdtdt  '•’Ojjoc^dtdt 


(26) 


When  the  above  double  Integrals  are  evalixated  for  two  different 
upper  limits,  say  t-^  and  three  sets  of  two  simultaneous 
equations  are  obtained.  Tne  constants  A,  B,  A',  B’,  Cp  ,  and 
Cdt  can  then  be  determined,  from  these  simultaneous  equations 
and  solutions  for  the  coefficients  Cm,  Cl  and  Cp  as  functions 
of  cx  will  have  been  foiind.  If  the  dynamic  pressure  q  can 
be  assumed  constant,  then  q  will  drop  out  of  equations  (24) 
through  (26),  and  the  Integrals  can  be  determined  by  numerically 
Integrating  the  cx  =  cx(t)  data.  If  q  Is  not  constant,  but 
quasi-steady,  then  q=q  (t)  may  be  determined  by  the  technique 
described  In  Method  VI,  and  the  above  Integrals  again  found 
by  nvimerlcally  Integrating  the  appropriate  functions  of  q(t) 
and  cx(t). 

l8.  Method  IV:  As  pointed  out  In  Method  III,  a  reasonable  assump- 
tlon  for  Intermediate  angles  of  attack  cx  Is  that  Cm=Acx  +  Bcx^ 
Substituting  this  expression  Into  equation  (3)  gives 


/ 


i 
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(Acx  +  A«Dm  =  lof 

A  manipulation  of  o<  yields 

(co»=^ cjCAc*  *bJ) 

where  cx)=^^  * 

Now  assume  q=  constant. 


If  Of  =  Ofj  and  a)  =0  when  t  =  0,  then 

This  equation  can  now  be  evaluated  for  two  different  upper 
limits,  say  cx.  j  and  Of,  >  giving  two  simultaneous  equations  from 
which  the  constants  A  and  B  and  hence  Cj^j,  can  be  determined. 

The  resulting  expressions  for  A  and  B  are 


A 

p-  I  a):fo«,-CTD-u>:^;-exu 

(«x;-cxp(of,-c<;  -  (pi;-ot‘)(o<t-o(i) 

The  quantltltes  a),,  and  cJ,  can  be  determined  by  fitting  a  curve 
to  the  C3(  =  o<(t)  data  using  any  of  various  methods,  such  as  the 

method  of  least  squares,  and  taking  the  first  derivative  of  the 
curve  with  respect  to  t.  This  technique  for  determlng  Cjvi  Is 
acciirate  only  If  the  dynamic  pressure  q  Is  essentially  constant. 
The  dynamic  pressure  can  be  determined  either  by  Method  VT  or 
Method  VTI. 


U);  (ext 


AmOr^q 
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19.  Method  V:  According  to  equation  (3), 


I  d 


r  -  ^  -  I  d  f,y\ 

qA.D„^‘PqA.Dl„  ^ 

Assxome  q=  constant.  Then,  if  Of  =0<i  and  ca)  =  0  when  t  =  0, 

]c>>do<  ^  gq/LO, 

CXj  O 

Now  let  cO  =0)0  when  0(  =  0.  Then 

fc^da  =-2^Ajb„ 

Since  two  models  are  used  in  each  shock  tunnel  test,  two  values 
for  the  term  on  the  right-hand  side  of  the  above  equation  can 
be  obtained  from  the  data  extracted  from  each  test,  if  both 
test  models  rotated  through  zero  angle  of  attack.  The  quantity 
oJo  is  determined  by  fitting  a  curve  to  the  Of  =  cx(t)  data 
around  the  region  where  0<  =  0  and  taking  the  first  derivative 
of  this  curve  with  respect  to  t.  Thus,  after  several  tests  are 
made,  several  values  of  the  term  -  lu)l/2qA„D^  can  be  determined, 
If  these  values  are  plotted  against  the  initial  angles  eXj in 
each  test,  then  a  curve  can  be  fitted  to  this  plot.  The  first 
derivative  of  this  curve  with  respect  to  cXjWill  then  give 
as  a  function  of  ofj,  or  angle  of  attack,  since 


r  - 

doTiV  2{^AJ 


k) 


Again,  as  in  Method  IV,  the  accuracy  of  the  above  method  of 
analysis  is  dependent  upon  the  requirement  that  the  dynamic 
pressure  q  be  essentially  constant. 

20.  Method  VI;  By  equation  (4), 

n  -JUsJb. 

Assume  that  q  is  quasi-steady,  so  that  for  relatively  short 
intervals  of  time,  q,  and  hence  X5,  can  be  considered  constant, 
then  fit  equations  of  the  form  Xl=C«s+C„t/  + Cesl/^  4o  short 

Intervals  of  the  Xg  =  Xg  (t)  data  by  the  method  of  least  squares. 
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If  Xg  is  nearly  constant  for  a  given  interval  of  time,  the 
second  derivative  x^  can  be  assumed  to  equal  Xg,  or  Xg=2C2s> 
and  the  dynamic  pressxire  q  In  this  Interval  can  be  accxirately 
estimated. 

21.  Method  VTI;  Assume  q  =  constant,  also  assume  the  velocity 
V,  of  the  sphere  Is  zero  and  Xs  =  Xoj^  when  t=0.  Then,  If 
equation  (4)  Is  integrated  twice  with  respect  to  t,  the  resulting 
equation  Is 


2*nsfXs-Xal) 


and  thus,  q  can  be  determined  directly  from  the  Xs  =  Xg(t)  data. 


22.  Each  of  the  above  methods  of  analysis  when  applied  to  the 
reduction  of  experimental  data  will  yield  results  of  varying 
accuracy  and  reliability,  depending  upon  the  type,  quality,  and 
quantity  of  the  data.  If  for  example,  exit)  Is  nearly  constant, 
then  Methods  I  and  II  will  give  the  most  accurate  results, 
whereas.  If  o<(t)  varies  considerably,  then  Methods  III,  IV,  and 
V  will  give  the  most  accurate  results.  Method  V  can  be  used 
successfully  only  If  several  tests  are  made.  By  use  of  Methods  I 
and  II,  a  single  value  for  Cj-,  C^,  or  Cm  Is  reduced  from  a  given 
Interval  of  the  data,  regardless  of  whether  or  not  the  angle  of 
attack  varies  during  this  interval.  Since  Cl  and  Cm  are  functions 
of  the  angle  cx  ,  an  average  value  oSf  is  determined  for  this 
Interval  of  data  such  that 


i“l 

where  Oi=  cx:(t),  and  n  is  the  number  of  data  points  In  the 
given  Interval.  When  a  single  value  for  either  Cl  or  Cm  Is 
determined  from  an  Interval  of  data,  the  quantity  5?  for  this 
same  Interval  Is  also  given.  For  models  that  are  symmetric  In 
the  plane  of  the  film,  the  drag  coefficient  Cp  Is  a  function  of 
the  absolute  value  of  the  angle  of  attack,  so  that  when  a  single 
value  for  Cp  Is  obtained  from  a  given  Interval  of  the  data  the 
quantity  for  this  same  interval  Is  also  given,  where 
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23.  A  major  source  of  error  in  the  results  stems  from  the 
difficulty  of  determining  with  precision  the  angle  cx  .  Though 
the  variation  in  o<  can  be  observed  with  accuracy  (roughly 

+  0.2°),  the  value  of  the  angle  of  attack  Itself  cannot  be 
^Tetermined  this  well.  This  is  due  to  the  effects  of  the  radial 
flow  in  the  4-ln.  Hypersonic  Shock  Tunnel  No.  3.  Though  the 
angles  of  attack  are  corrected  for  flow  radlallty,  they  are 
corrected  only  to  the  extent  of  the  accuracy  to  which  the  degree 
of  radlallty  is  known.  The  accuracy  of  the  measurement  of 
is  best  at  the  high  Mach  nxambers,  since  the  radlallty  of  the  flow 
is  least  at  these  positions.  This  error  in  ex  will  introduce 
greater  error  in  the  results,  if  the  results  are  obtained  by 
numerical  Integrations  of  the  cx(t)  data,  such  as  is  the  case 
if  Method  III  is  used. 

24.  Ordinarily,  for  a  series  of  several  tests,  the  scatter  of  the 
results  can  be  kept  within  +  5  percent  of  a  mean-value  curve 
(except  for  the  very  small  values  of  Cl  and  Civj  such  as  those  that 
occur  at  small  angles  of  attack),  but  if  the  test  models  are 
relatively  massive  (e.g.,  greater  than  15  grams  for  a  model  3 
Inches  in  length),  then  their  accelerations  will  be  small,  and 
the  given  limit  of  scatter  in  the  results  will  be  surpassed. 

25.  It  is  the  opinion  of  the  authors  that  a  model  should  be 
tested  at  least  twice  at  roughly  equal  conditions.  It  is  felt 
that  the  Increase  in  reliability  of  the  results  will  more  than 
offset  the  cost  of  additional  tests. 
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APPENDIX 


THE  METHOD  OP  LEAST  SQUARES 

1.  Given  Is  a  set  of  data  (xj^,  t^)  for  which  x  =  x(t)  and 
i=l,2,3> . . • jn.  It  is  desired  that  a  curve  of  the  form  x'=Co 
+Cit+C2t^  be  fitted  through  these  data  in  such  a  way  that  the 
sm  of  the  squares  of  the  individual  deviations  x^-x'  be  a 
minimum,  or 


G=X[xi-(Q.+  Cjt,+ C,t:)] 

1=1 

be  a  minimum.  Now,  G=G(Co,  C2)>  so  that 

dG=ifdc..|ed,^.|f^dc. 

Then  the  conditions  for  the  fitted  curve  will  be  satisfied  when 

ic ) 

“  i=l 

ic,  (0.+ c,t,  +  c,t',)T=o 

ic,  (g.*  <-,1,+  C,t')J=0 

Upon  carrying  through  the  above  differentiations,  we  have 


-2^JXi-(C„+C,ti+  c,tj)]  =  o 

i'l 

n 

-2^t,[Xj-  (&+C,t.+C,t-)]=0 

1=1 

i=l 
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or 


Zx,=  C.n+c2^ti  +cXt! 

i=l  i=l  i=l 

n  n.  n.  n 

Z  Xitr  cXti 

i=l  i=J  i=l  i-1 

+cXt“ 


oe 


The  problem  Is  now  one  of  solving  these  three  simultaneous 
equations.  This  Is  done  most  simply  with  the  use  of  determinants. 


2.  Let 


D  = 


n 

n 

Zt. 

n 

Zt: 

i-1 

n 

n 

ztr 

i=i 

n 

Zt: 

H 

i=I 

n 

ztr 

i^l 

n 

Zt: 

1=1 

i-  1 

n 

Zt, 

i=i 

n 

Zt: 

i=l 

n  n 

Ex,t.  ztr 

:i=i  1 1 

n 

Zt: 

i--i 

n 

Zt’ 

i=i 

n 

Zt: 

i=l 

n 

n 

Zx. 

i=l 

n 

Zt: 

1=1 

n 

Zt. 

i=l 

Zx.t. 

ri  i-1 

n 

Zt! 

i=l 

n  n 

Zx.t.^  Ztl 

i=l 
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ri  n 


n 

Zt, 

in 

Zx. 

i=l 

n 

>;t, 

1=1 

n 

Etf 

i-l 

i-1 

n 

ztr 

rl 

n 

i=l 

Then  finally. 


-a- 

D 


"  D 


3.  If  It  Is  desired  that  an  equation  of  the  form  x'=Co+ 
be  fitted  to  the  same  data  by  the  same  method,  then 


n 


c,tj' 
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and  by  the  same  treatment. 


c= 


Xt. 

i  x.t,  z  t: 


n 


Zt, 

id 

n  n 

zt  ztr 


Ca= 


n 


n 


n 

n 


n 

ZXi 

i=l 
n 


Zt, 


n 

Zt. 

i=l 

Zt.  Zt; 

id  i=l 
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TEST  SECTION 
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FIG.  3  SCHEMATIC  DIAGRAM  SHOWING  LOCATIONS  OF 
WINDOWS  IN  RECEIVER  TANK  OF  N.O.L. 

4  inch  HYPERSONIC  SHOCK  TUNNEL  NO.  3 
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FIG.  4  VIEW  OF  MODELS  FROM  INSIDE 
RECEIVER  TANK. 


CONDITIONS  IN  THE  NOL  4 -IN.  HYPERSONIC  SHOCK  TUNNEL  NO. 
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